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Abstract 

Let A be the limiting smallest eigenvalue in the general (/3, a)-Laguerre ensem- 
ble of random matrix theory. That is, A is the n f oo distributional limit of the 
(scaled) minimal point drawn from the density proportional to rii<j<j<n 1^* ~~ *M 

nr=i\^ (a+1)_le "^ Al on ( R +)™- Here /3 > 0, a > -1; for (3 = 1,2,4 and integer a, this 
object governs the singular values of certain rank n Gaussian matrices. We prove that 

8 I — t("/+1 — /3 /2) 

P(A > A) = e-f A + 2 ^A ~P e(/3, a)(l + o(l)) 



as A t oo in which 7 = ^(a + 1) — 1 and e(/3, a) > is a constant (which we do not 
determine). This estimate complements/extends various results previously available 
for special values of (3 and a. 



1 Introduction 

The shape of the distribution of the smallest singular value of a "typical" matrix is a deeply 
studied question. An overview of the varying motivations for this problem may be found 
in [TT]. In the case of Gaussian matrices, many exact formulas are available both at finite 
dimension and asymptotically 01 [T3] . Only quite recently has it been shown that the asymp- 
totic laws are universal beyond the Gaussian case (in the sense of being insensitive to the 
statistics of the matrix entries), see |12j . 

Here we consider the "general beta" analogues of the classical Gaussian ensembles. These 
are defined by placing a measure on n nonnegative real points Ai, A 2 , . . . , A n with density 
function (a normalization constant times) 

n lA.-A/nAf" 11 - 1 ^-. (i-D 

l<i<j<n i=l 
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When (3 = 1, 2, 4 and a = 0, 1, 2, . . . this is the joint square-singular value law of an nx (n+a) 
real, complex, or quaternion Gaussian matrix. It is however a sensible law for any /3 > 
and a > — 1, and, what is more, still a joint square-singular value law for a certain random 
bi-diagonal matrix ensemble [3]. Further, the least order statistic \ m in satisfies a limit law: 
as n I" oo, n 2 X m i n converges (in distribution) to a well defined random variable, denoted here 
by A (= A(/3, a)). There are several proofs of this for special values of /3 and a; [10J contains 
a proof (making use of the bi-diagonal representation of [3] and substantiating a conjecture 
of [5]) valid for all values of those parameters. 

Our starting point is a relation between the law of A and the explosion/non-explosion of 
the diffusion process: with b a Brownian motion, 

dx(t) = db(t) + (^{a + i) - ^V\e~P t/s cosh x(t)^J dt. (1.2) 

In particular, a corrected version of Theorem 2 of [10] (see also the derivation leading to 
(12. 3p below) implies that 

P(A > A) = Poo, (t H> ar(t) never hits - oo). (1.3) 

Here P CjS indicates the law on paths induced by x, begun from position c at time s. Our 
main result reads: 

Theorem 1. Let p\ = pA,/?, a denote the right hand side of U.3\) . For large values of A it 
holds 

p x = e -fA+27>/A; r 7(T V /a) e (/3, a )(i + o(l)). (1.4) 
Here 7 = |(a + 1) — 1 and e(/3, a) > is an undetermined constant. 

There has already been a great deal of work in this direction, though focussed on dealing 
directly with the statistics (11. ip rather than our passage time description (II. 3p . The funda- 
mental treatment of Tracy- Widom [13] for /3 = 2 produced the correct A — > 00 asymptotics 
of pA.2,a up to a multiplicative constant and provided a conjecture for that constant, e(2, a). 
This has recently been verified by Ehrhardt [7], for \a\ < 1, by operator theoretic techniques, 
and for all a > —1 by Deift-Krasovsky-Vasilevska [2] using Riemann Hilbert Problem ma- 
chinery. A non-rigorous argument in [1] predicted all factors in the asymptotics save the 
constant for all (/3,a). Making use of integral identities available at special values of /3 and 
integer a, Forrester has a sound conjecture for the value of the general constant e(/3,a), see 
[6] . It appears possible that a type of analytic continuation argument could extend the result 
of [6] to all (3 (though for still integer a). The method employed here leaves c(/3, a) in opaque 
form, as a somewhat involved expectation over diffusion paths; an explicit determination of 
this object for all and a remains an open problem. 



2 



In many ways, the chief insight of this paper is to cast the diffusion (11. 2p . which encodes 
the desired probability distribution, in the present form. (The process which appears in 
[TO] is related by a change of variables.) In fact, t i— > x(t) is remarkably similar to the 
process studied by Valko-Virag in estimating the probability of large gaps in the general 
beta "bulk" [T3]. They showed that the probability of a gap being larger than A is equal to 
the non-explosion, again to — oo, of 

dz(t) = db(t) + Qtanh^(t) - ^Ae _/3t/4 cosh z{t)j dt, 

begun again at +00. It is no surprise then that their basic argument, which involves esti- 
mating the Cameron-Martin-Girsanov factor produced by a well-chosen change of measure, 
may be followed in this case. 

The proof of Theorem 1 occupies sections 3 and 4; section 2 gives a self-contained expla- 
nation of the identity (II. 3p . 



2 Passage time description for A 

Without pointing the reader to [10] and the subsequent erratum, it is easy enough to give a 
brief derivation of the relevance of the diffusion (II. 2ft to the distribution function P(A > A). 
The main result of [1QJ shows that A -1 is the maximal eigenvalue of the almost surely trace 
class integral operator 

poo / ptAs \ 

Lf,, a il>(t) := J (J e au+ ^ b(u) duj ^(s)e- {a+1)s -7? b{s) ds, (2.1) 

acting on L 2 [R +J p], /i(dt) = e ~ (a+1)t ~^ m dt. Here t \-> b(t) is a standard Brownian motion. 

Any nonnegative L 2 solution of ip{t) = \Lp ja iJ;(t) satisfies if)(0) = and ip'(t) : = 
dif)(t) /dt > for all t > 0, as can be seen by taking derivatives of both sides of the eigenvalue 
equation: 

^'(t) = \ e at+ VP b ^ J™ tfj( s )e- {a+1)s -VP b{s) ds. 
This converts to a differential system: 

dijj'{t) = -=ij/(t)db(t) + [{a + -z)ij/{t) - Ae _ V(t)]*, = ^'{t)dt (2.2) 

VP P 

which can be used to test whether a fixed A is at or below an eigenvalue. Specifically, A 
is strictly below the groundstate eigenvalue A if the solution to (12. 2p begun at tp(0) = 
(and ^'(O) = 1 say) satisfies ip(t) > 0,ip'(t) > for all time (note that solutions of (I2.2p are 
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decreasing in A). It is now the standard trick to translate this condition onto the diffusion 
q(t) := i/)'(t)/ij)(t) which solves 

dq(t) = -^=q(t)db{t) + {{a + l)g(t) - q\t) - Xe^dt, 

started from +00 at time t = 0. In particular, if r c is the passage time of q to a level c, the 
event {A > A} coincides with {r = 00}. Now the change of variables, 

x(t) := log(?(#/4)) + 0t/8 - log A/2, (2.3) 

explains the identity (II. 3ft . 

As a bit of amplification, we remark that for q = q(-; a, (3, A) with a > 0, 

P(r_ 0O (g) < oo|r (g) < 00) = 1. (2.4) 

So, at least for a > 0, one can replace the condition of q never vanishing with the (more 
familiar) condition that q never explodes to —00. Furthermore, a change of variables similar 
to (12.31) shows that the event that q(-\ a, (3, A) started from never hits —00 is the same as 
the event that q(-; —a — I, (3, A), started from +00 never hits 0. In particular, defining instead 
x(t) : = -log(-g(/3t/4)) + f3t/8 - log A/2 (keeping in mid that for all t > r , q(t) < 0), x(t) 
will solve (11.21) with a replaced by —a — 1. One concludes that 

lim P(A > A) = lim Poo(r (g(-; a, (3, A)) = 00) = limP (r_ oo (g(-; a, 0, A)) = 00) = 

aj.— 1 ai— 1 afO 

for any A > (by say monotone convergence), as would have been guessed ahead of time. 

To prove 02.4)) . on the event {r < 00} introduce the simpler change of variables u(t) = 
log(—q(t + r )). This process satisfies 

du(t) = -^=db{t) + {a + e u{t) + Xe'^e-'e-^dt, u(0+) G (-00, 00), 
VP 

to which we compare the homogeneous process defined by 

dv(t) = -^=db(t) + [a + e v(t) }dt, v{0) = u(0+) G (-00, 00). 

VP 

As w(t) > f(t), q explodes to —00 in finite time if v explodes to +00 in finite time (we 
continue to work on the event {ro(q) < 00}). 

Now apply Feller's test, in the form given by Proposition 5.32 (part (ii)) of [9j. In 
particular, bring in the Lyapunov function 

[ x f y 1 {Pf x 

mix) = / s(y) / ——dzdy where six) = exp / (a + e z )dz 

Jo Jo s{z) V. 2 J 



(s(x) is the derivative of the scale function for v). Since 



lim m(x) < oo, while, if a > 0, lim 771(2) = —00, 

X— >CO X— >— OO 

the cited form of Feller's test implies that S — f{t : v(t) ^ (—00, 00)} is finite with proba- 
bility one. However, it is impossible that v(t) ever hits —00 (it is easily bounded below by 
a Brownian motion with constant drift a). This completes the proof. 

3 Change of measure 

Hereafter it is convenient to put the time index in subscripts, i.e., x(t) becomes Xt and the 
like. To begin, introduce the notation 

px(c) = F c (x t never explodes). 

Then, by the strong Markov property, 

p A = p A (oo) = pi(x T ),Xi > -00 for t G (0,T] , (3.1) 

upon choosing 

T = |logA. (3.2) 
The change of measure is now enacted on the expectation (13. ip . 

Proposition 2. Let h(t, x) be C 1 in both variables and bounded for t < T. Then, the law 
on paths up to time T induced by 

dy t = db t + (h(t, y t ) - ^e'^ 8 sinh(7/^ dt, y = 00 

is absolutely continuous with respect to that oft^ x t , xq = 00, subject to x t > —00, < t < 
T. Moreover, 

pA = E oc [p 1 (7/ T ) J R T (7/.)], (3.3) 

in which, for s < T, 

\ogR s (y.) = J (f(t,y t )-g(t,y t ))dy t -^J (f(t,y t ) - g 2 {t,y t ))dt, 

f(t,y) = f (a + |) - fVAe-W 8 coshy and g(t,y) = h(t,y) - |v^e _/3t/8 sinhj/. 

This is just the formula of Cameron-Martin-Girsanov, applied to the particular case of a 
diffusion with explosion for which it is important to point out that the test function Pi(xt) in 
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question vanishes when T is larger than the explosion time. One also notes that the general 
form of the y-drift, effectively a bounded function minus sinhy, allows y% to be started at 
+00 and prevents y t from exploding on [0, T]. To then carry out the standard proof of 
Cameron-Martin-Girsanov in the present context, it must be checked that R^~ 1 (x), which is 
a local martingale by construction, is actually a martingale. But again by the general form 
of the y-drift, both f — g and f 2 — g 2 are bounded when the path is bounded below, keeping 
R^ 1 (x) bounded prior to the explosion time of Xt- Plainly, R^~ (x) = at and after the 
explosion time. 

The first ingredient of the proof of Theorem [1] is the following. Throughout the below, 
[x]~ denotes the negative part of x 6 R. Also recall that 7 = -|(a + 1) — 1. 

Lemma 3. There exists a choice of h in Proposition 2 so that, for appropriate u, <ft satisfying 
W(t, y)\ < Ki + H2[y\~ f or allt > and constant K\, k 2 , and \4>(t, y)\ < 4>{t) with f °° <f)(t)dt < 
00, it holds that 



logJfrfo.) = -f A + 2 T VA- ( 7(7+ 2 1 ~ /3/2) ) logA 

+ v(T,yr)+ [ <j>(T-t,y t )dt. 
Jo 



(3.4) 



+ -e~ v < 
2 



Once h is in hand, the lemma is readily verified. In particular, 



h(t, y) = |(a + I) + hi(y) + e-^h 2 {y) (3.5) 



where 



HV) = - T ^, (3-6) 

hiiy) = — !^ {ih 2 M - h 2 M) + § (a + \){hM - h(0)) + (h[(y) - h[(0)) 

That both hi and hi are uniformly bounded, hi being integrable at +00 while h 2 is integrable 
at both ±00 figure into the bounds on v and </> in the lemma. 

It is more instructive however to describe how h is discovered, each step achieving suc- 
cessive order in A, vA, etc., and the various bounds claimed in the lemma seen along the 
way. 

Step 1 begins by expanding out the exponential R? factor with a generic h: 

32 pT 1 pT o2 



\ogR T (y.) = -L X ^ e-Mdt + lJ tf(t,y t )dt-^a+ l -) 2 T (3.7) 

e-^ 8 h(t,y t )smh(y t )dt + ^-(a + h^/\ f e~^ 8 cosh(y t )dt 

8 1 Jo 



2 

f VX I* e-W 8 e-y<dy t - f [h{t,y t ) - |(a + \)]dy t . 
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By the choice of T, the first term equals — f (A — 1) which already gives the leading order 
and explains the particulars of the sinh?/ term in the choice of the ?/-drift . The last term, 
coupled with the fact that y = oo, prompts a natural shift of h by the factor f (a + |). That 
is, h is replaced with h + f (a + \). 

Step 2 enacts the above shift, and also introduces the obvious ltd substitution in the 
second last term of (13.71) . 

2 Jo 2 4 16 Jo 



to write: 



?2 rT 



log R T (y.) = — — A / e'^dt (3.8) 



+ (~ £e-W*e-*dt - J e-W 8 h(t,y t ) sinh(^) 



+ h 2 (t,y t )dt+^(a+^) J h(t,yt)dt-J h(t,y t )dy t + ^e~ yT . 



This draws attention to line two of (13.81) . which should produce the final constant times vA 
term. This may be achieved most easily by introducing a deterministic integrand in that 
line via the substitution 

w*, y) = l e —^ + kfa y) ■= h ^y) + n ^ v)> ( 3 - 9 ) 

2 smh(y) 



so that 



T r T 



J / e-W 8 e- yt dt- / e-W 8 h{t,y t )smh{y t )dt 



T rT 

e -W*dt- / e-P t/8 h(t,y t )smh(y t )dt. 



o 



Evaluating all deterministic factors thus far, step 2 is summarized by 

log R T (y.) = -|a + 2 7 v/A + t e -vr - (/3(a + I) + 2) (3.10) 



~V\f e-W 8 h(t,y t )smh(y t )dt 



+ 2J h (^yt) dt +j( a +2^J h (^yt)dt-J^ h(t,yt)dy t . 

The first two terms above exhibit the proposed order A and order a/A factors in the statement 
of the lemma, showing that there was not much flexibility in the choice of the (uniformly 
bounded) function hi in (13. 9L 
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Step 3 is to pin down the log A factor in the exponent (or, equivalently, the T factor). A 
look at line two of (13.101) suggests a prescription for h: 



h(t,y) 



-G8/8)(T-t) 



h(y) 



(3.n; 



in which /i 3 must be chosen so that h 2 is bounded (and more). 

With r](t) = |e _/3 */ 8 , we employ Ito's lemma once more to write the final term in (13. lOj) 



as m 



h(t,y t )dy t = H 1 {y T ) + H 2 (T,y T ) 



(3.12) 



r 



h[(y t )dt 



T pyt 

r}'(T-t)[ h 2 (z)dz]dt 



rj(T-t)h' 2 (y t )dt. 



Here Hi and H 2 denote the anti-derivative terms which appear: Hi{i/t) = J * h\{y)dy L and 
H 2 (T,yT) = r)(T — t) h 2 (y)dy |^ . Note that the boundary values of H 2 will necessitate 
that our choice of h 2 , like that of hi, is integrable at +00 (= y ). Now expand out the last 
two lines of ( 13 . 10[) to read: 



t r 



\h\{yt) + \h\{y t ) + j{a + ^)hi{y t ) - h 3 (y t ) 



dt + Hi(y T ) + H 2 (T,y T ) 



(3.13) 



+ 



T 



-V 2 (T-t)ht(y t ) 

( 1 3 f yt 1 

+ n{T - t) I hi{y t )h 2 {y t ) + J (a + -)h 2 (y t ) + |[ / h 2 {y)dy\ + -h' 2 (y t ] 



dt. 



The first term of ( 13. 131) prompts the choice of namely set 

Ms/) = \ h l(v) + \ h 'M + j (a + ~)MsO - «, 

for a constant k which makes h 2 (y) = h 3 (y)/ sinh(y) bounded. We find that 

k = \hl(0) + \h'M + I (a + i)M0) = ^^f^ - f (a + 5)7 = -|(7 + 1 - P/% 



compare (I3.6p . In other words, with this choice the first term of ( 13 . 1 3[) equals ^/clogA, the 
advertised log A contribution in Theorem [TJ 



To finish the proof of the lemma we further identify 

u(T, y) = -(P(a + l -) + 2) + H^y) + H 2 (T, y), 



(3.14) 



and <p(T — t, y) with the integrand of the last term in (13.131) . One now checks: hi and h 2 
along with their derivatives are uniformly bounded over the entire real line (with constants 
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depending on a and (3 of course), hi is integrable at +00, and h 2 is in fact integrable at 
both ±00. The boundedness and integrability (at +00) of h\ and h 2 give immediately that 
l-^i (v)\ + \H2{t,y)\ < K\ + ^2U -i proving the claimed bound on v. For <ft one shows that 
\4>(T— t, y)\ < cr)(T—t), using the additional appraisals on the derivatives, that \h 2 \ < 00, 
and the simple fact < i](T — t)< i](0) for t 6 [0, T}. 

4 Constant term 

The conclusion of the previous section is that 
with 



p 1 ( 2/T )e§ e ' !/T + l '^)+/o T ^-**H ; 

and v and satisfying the bounds outlined in Lemma [3) It remains to show that the existence 
of a (non-zero) constant e = c(a, (3) such that lim^oo t x = e. This is again structurally 
identical to [14J. 

The first observation is that the Eoo integration is performed over paths that are mono- 
tonically decreasing in T. The nicest way to see this is to replace the integration over 
Dti < t < T with that over 

yf = y t+ T, -T<t<0 

which satisfies 

dyj = dh t + (h(t + T, yj) - ^e~* sinh yj)dt, y T _ T = 00. 

If this family of processes is run on the same Brownian motion, t 1— > b t , it follows that 
yj 1 < yf 2 f° r t > by definition y^} T < y^ T and the evolution maintains the ordering. 

Denote this sequence of corresponding expectations simply by E and record that 

e A = E[ Pl (y T )e^ T >], iP(y T ) = ^ + „(T,j£) + J%{t,y T _ t )dt. (4.2) 

Next, pick a constant h$ such that 

inf hit + T,y)> h 

-oo<j/<oo,-T<t<0 

(a look at (13. 5 p and (13 .6p shows this is possible), and introduce the stationary diffusion 1 1->- z t 
on the negative half-line with generator 

1 d 2 d (3 

£=-——+ f(z) — , f(z) = h sinhz, 

2dz 2 JK ' dz JK ' 2 



9 



and reflected (downward) at the origin. In particular, for all t > —T, F(z t G dz) = m(dz) 
where 

m(dz) = KQ e 2h ° z -P coshz dz, z G (-00, 0], (4.3) 

and Kq is the appropriate normalizer. This is the well-known formula for the speed measure 
(see for example [5]), or one may check that J"° C(j)(z)m(dz) = for all smooth satisfying 
0'(O) = 0. 

Again running z t on the same Brownian motion, it holds that yj > zt > —00 for all 
t G [-T, 0]. This is plain at the starting time, and continues by the domination (from below) 
of the y T -drift by that of z. It follows that there exists a random variable y^° > —00 such 
that 

lim yj = y^° almost surely for each t G (—00, 0]. (4.4) 

T— >oo 

To pass this convergence under the E-expectation we prepare the following (and defer the 
proof to the end of the section). 

Lemma 4. The function x pi(x) is continuous, strictly positive on x > —00 and satisfies 

p 1 (x) < K 3 e~* e ~ x (4.5) 

for a constant k 3 . 

Courtesy (14.41) and the first statement in Lemma H] we have that 

lim Pxiy^e^ = pMe §e-^ + u(oo,y^ + f-^t,y-) d t ( 4 . 6 ) 

:= Pl ( 2/0 oo )e^^), 

using continuity (for the first three factors) and dominated convergence (for the last factor). 
The evaluation i/(t, y)\ t=00 simply has the effect of setting of one if2-terms of which v is 
comprised to zero, recall f )3.14p . 

Next, by the estimates on u, <fi from Lemma [3] and (14. 5p . there are the bounds 

^(y^e-^r < Pl (yT)e^ < n A e^' + ^ , (4.7) 

with positive constants k 4 , k 5 . 

Note that both bounds in (14.71) depend only on the marginal of the process at time 0, 
and denote the left and right hand sides by p_(y^) and p + (yj) respectively. Invoking again 
the path-wise control, yf > z t , t G [— T, 0] we have that 

Pi(yo)e^ iyT) < p+(*o), E[p + (^ )] = I P+ (z)m(dz) < 00, 

^—00 
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where m is defined in (14. 3p . Returning to (I4.ip . ( 14. 6 j) and dominated convergence now 
produce 

lim e A = lim E[ Pl (y T )e^ T )] = E[p 1 (^)e^^°°)] := e, 

A— >oo T->oo 

defining the constant e in the statement of Theorem 1. That c is not equal to zero follows 
from 

e > liminf Efpi^e^ >] > / p_(z)m(dz) > 0. 

T ^°° J-oo 

Here we have used that z !->■ p_(^) is decreasing in order to replace y T -paths with z-paths, 
along with the fact that p\(z) (and so too p-(z)) is strictly positive (Lemma 0J. This 
completes the proof of Theorem 1, granted the below. 

Proof of Lemma 4- The continuity follows from that of the transition density p(-,x,y) in 
both space variables (the corresponding generator is hypo-elliptic). 

To see that pi(x) > 0, first note that the operator Lg iQ defined in (12 .li which encodes 
the point process of eigenvalues is positive and compact. A proof that Lg )0 is in fact (almost 
surely) trace class is contained in Lemma 6 of [10J. Its maximal eigenvalue, A -1 , is therefore 
almost surely bounded above, and so there exists a small enough A > such that p\ = 
pA (oo) > 0. Next, by the Markov property, 



poo 

Pa (oo)= / p{t,oo,x)p Xoe -0 t /4(x)dx, 



and it follows that for every t > there is a x such that p\ oe -pt/4.(xo) > 0. Applying the 
same formula once again, we find that for any zeK 

/oo 
P(s,x,y)p e -0s/i(y)dy > F x (x s > z)p e -p./*(z). 

To finish, choose s = t — 4 log Ao and then set z to be the appropriate xq. 

For the bound (14. 5 p we can restrict to x less than some large negative constant, and note 
that pi(x) is bounded by the probability of non explosion for the following process 



dy t = db t + ^a + ^-^e-W 8 e-*. 



since the downward drift on y is weaker than that of x. Next make the change yt = yt + y 
to obtain the homogenous process 

dyt = db t + ^(a + l- e- yt ) , 

to which we can apply a version of Feller's test, similar to what was done at the end of 
Section 2. A scale function for the ^/-process is 



s(v) 



J\xp{-^[(a+l)^ + e^-l]}d^ 
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and the probability that this process exits through +00 is exactly the probability of not 
exploding. This follows by checking the conditions of now Proposition 5.22 of [9J. According 
to that same proposition, the exit probability equals 



six) — s — 00 



IjT aq>{-f [(a + l)e + e-^]}de, 



s(+oo) — s(— 00) Z 

from which the required bound easily follows. □ 
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